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The field radiated during the current surge* at the source is derived by a new method. 
The; effects of changing such parameters as rate of current rise and speed of channel develop- 
ment are indicated. Amplitude spectra at 100 km an; deduced for typical first and subse- 
quent strokes in a discharge carrying negative electricity to earth, and for a positive stroke. 
The three spectra peak at approximately 5, 3, and 3 kHz, respectively. 

Mean amplitude spectra at 100 km are also deduced by extrapolation back, making 
allowances for propagation, of atmospherics observations. The average spectrum for a 
group of atmospherics obeys a simply normalized relation involving the frequency of maxi- 
mum amplitude; this frequency changes between groups over the range of 4 to 7 kHz. 
The amplitudes, both on a broadband and a single-frequency basis, are log-normally 
distributed, with standard deviations of some 7 dB. 

Averaged spectra for groups of observations agree quite well with those computed 
from the stroke parameters. However, it is found that the spectral shape varies extensively 
among individual atmospherics. 



1. Introduction 

A companion paper [Arnold and Pierce, 1964] lias 
surveyed the generation of VLF atmospherics during 
lightning flashes and has considered in detail the 
amplitude spectra associated with the L pulses accom- 
panying stepped leaders and with the K changes 
occurring during certain stages of both cloud and 
ground discharges. The present paper examines the 
VLF atmospherics radiated during the return stroke 
of a flash to earth. These atmospherics represent 
the largest VLF disturbances associated with light- 
ning but, as indicated in the companion paper, are 
not the sole VLF signals of importance. To a first 
approximation, a return stroke can be considered as 
a discharge in a vertical channel, growing upward 
from the ground to neutralize the charge deposited 
by the preceding leader. If the form of the current 
surge and the rate of growth of the channel are 
known, the spectrum of the resulting atmospheric 
can be estimated by applying Fourier methods to 
the second time derivative of the electric moment 
[Lejay, 1926]. The electric moment M is usually 
defined as twice the summed product of the charges 
and their heights above ground, the factor two taking 
account of the image charges. The "current mo- 
ment" is often used by radio engineers; it involves 
the product of current and the length over which 
the current flows, and is approximately equal to 
dM/dt 



2. Theoretical Approach 

2.1. Background 

As the treatment of return strokes by Bruce 
and Golde [1941] summarized rather well the observa- 
tions made up to that time, and as their model 
for a return stroke has been widely accepted, a 
brief description is desirable before modifications, in 
the light of more recent work, are discussed. Bruce 
and Golde considered the current surge at the ground 
to be of the double-exponential form 



i t =io[exp (—at) — exp (—#)], 



(2.1) 



where i t is the current at the base of the return 
stroke at time t, and i , a, and are constants. 
They suggested as typical values 30,000 A, 4.4X10 4 
sec" 1 , and 4.6X10 5 sec -1 for i , a, and f3, respectively. 
From photographic data, notably that of Schonlaiid 
[1956], Bruce and Golde concluded that the first 
return stroke of a flash to ground grows upward 
at a rate given by 



y,= y ex P (-yt) 



(2.2) 



where V is 8X 10 7 m sec' 1 and y is 3 X 10 4 sec" 1 . For 
subsequent strokes in a lightning (lash, the upward 
growth is relatively constant, so that y is equal to 
[Schonland. 1956]/ 
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Bruce and Golde give a simple treatment of their 
return-stroke model in which it is assumed that the 
current in the channel is uniform. This requires 
the instantaneous transfer of charge from the base 
of the stroke to the tip. Hence 



dM r< 

dt 2%t J 



V t dt- 



2i V c 



[exp (—at) 



-exp (-#)][i- eX p (-yt)}, (2.3) 



where M is the dipole moment. (The factor 2 allows 
for image charges.) [Lejay, 1926.] Equation (2.3) 
can be differentiated with respect to time and then 
treated by Fourier methods to obtain the source 
spectrum for both first and subsequent strokes. 
The resulting spectra peak near 9 kHz for first 
strokes and 7 kHz for subsequent strokes [Pierce, 
I960]. 

2.2. A Modified Treatment 

A wide range of spectra can be generated by 
varying the Bruce-Golde parameters for return 
strokes. Furthermore, additional degrees of freedom 
can be introduced by abandoning the assumption 
of uniform current in the channel, by permitting 
variations in the charge laid down per unit length 
by the leader stroke, and so on. The additional 
degrees of freedom greatly complicate the derivation 
of the source spectrum. As the spectrum of an 
individual stroke is influenced by such things as 
nonvertical sections of channel, branching, and 
irregularities of the current surge, which cannot 
be satisfactorily included in any simple return- 
stroke model, a point is soon reached beyond which 
further elaboration of models is of little practical 
value. 

The following treatment avoids the unrealistic 
assumption of a uniform current and, hence, instan- 
taneous charge transfer, in the lightning channel. 
Consider a first return stroke, described by (2.1) and 
(2.2), with the charge in the channel moving at a 
constant velocity u, where u> V . If h t is the chan- 
nel length at time t, obviously 



■/. 



h t = V t dt. 



(2.4) 



The charge arriving at the tip is "frozen" there to 
neutralize charge of opposite sign laid down by the 
leader stroke. Thus, at any time t, the channel 
will contain a flowing charge of density p per unit 
length and be surrounded by a frozen charge of 
density p' per unit length. The quantity p is a 
function of height, h, and time, t; it is given by 



{-K-SHhH)]}- 



(2.5) 



Let p t be the density of flowing charge at the tip. 
Then the charge laid down at the tip per unit time 



is p*(u— V t ), and the density laid down per unit 
length is 

v; 
Pt=Pt 



where 



ra 



^M-(«-$HK'-S)]}- 

(2.6) 
The electric moment at time t, M t , is given by 

M,=2 £' ( P + P ')h • dh, 



(2.7) 



where the factor 2 allows for image charges. By 
standard methods for differentiating definite inte- 
grals, and recognizing that p 1 does not depend on t y 



dM t 

dt ~ 



'■£ 



Jo ^ h ' dh+2 \ P-^7-dh+2(p t + p' t )h t ^ 



dt 



dt 

(2.8) 



As h is not a function of t, the second term on the 
right is zero. Differentiating (2.4) and (2.5) and 
substituting into (2.8) lead, after some manipulation, 
to 

dM t 2tyM , , N f fah\ .,"1 

-^exp(-0»[exp(<|)-l]. (2.9) 

Differentiating again with respect to time yields 

d 2 M . , ,, r MAI 

— =2 ^exp(-aO[l-exp^-Jj 

+2i y exp (- 7 exp [-a (*~)] 

-2i u exp (—#) Tl-exp (^)l 

-2i V exp (-7O exp [-£ (*~^)} (2-10) 

It is seen that, for u equal to <» , (2.9) becomes identi- 
cal with the Bruce-Golde expression for dM t /dt, as 
would be expected. 

2.3. Influence of the Various Parameters on Return- 
Stroke Spectra 

In order to study the effects of varying the param- 
eters in (2.10), two computer programs have been 
set up. The first provides d 2 M/dt 2 at intervals of 
1, 2, 3 . . . 100 /xsec; the second computes d 2 M/dt 2 
at intervals of 2, 4, 6 . . . 1000 /zsec, and then uses this 
output to derive the Fourier coefficients at frequen- 
cies of 1, 2, 3 ... 25 kHz, By inserting some 50 
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combinations of the parameters into the program, 
the effects of changes in the individual ones have been 
determined. Results have been normalized to show 
electric field components at 100 km from the source. 

Figures 1 and 2 show the effects of changing u 
from oo to a finite value. Curve (a) of figure 1 
shows d 2 M/dt 2 as a function of time, t, for a first 
return stroke with the following parameters: a = 
4.5X10 4 sec" 1 , 0-4.5X1O 5 sec" 1 , 7=3X10 4 sec" 1 , 
^ =8X10 7 m/sec" 1 , ^=30,000 A, and w= «>. Apart 
from very minor changes in a and 0, this is the Bruce- 
Golde first return stroke, with a maximum current 
near 20,000 A. Curves (b) and (c) of figure 1 show 
<PMIdt 2 with u reduced to 3X10 8 m/sec" 1 and 8X10 7 
m/sec -1 , respectively, with all other parameter sun- 
changed. The speeds chosen, that of light and of 
V , represent upper and lower limits for it. Curve 
(a) shows a positive peak near 3 jusec and the maxi- 
mum negative excursion at 35 /xsec. Little change 
is evident on curve (b), apart from a slight delay 
of the negative excursion. However, reduction of 
u to 8X10 7 m/sec -1 (curve (c)) moves the positive 
peak to 5 /xsec and the greatest negative excursion 
to 43 /xsec. Curves (a), (b), and (c) of figure 2 
show the corresponding spectra; the spectral peak 
(/'„,) moves from 9.5 kHz for u= °° to 9 kHz for u= 
3X10 S m/sec -1 , and to 8 kHz for ?x = 8X10 7 m/sec" 1 . 

The maximum length attained by a first stroke 
channel is given by V /y. The Bruce-Golde stroke 
has a maximum length of only 2.7 km, which seems 
low for some parts of the world. It would appear 
that lengthening the channel would lead to a lower 
value for f m . The lengthening has been done in the 
computer program by reducing 7 to 2X10 4 sec -1 , 
which leads to a total length of 4 km, and by using 
all other parameters from curve (b) of figure 1. 
Figure 3 shows d?M/dt 2 for the lengthened stroke as 
curve (a): Curve (b) (7=3 X10 4 sec -1 ) is reproduced 
from figure 1 for comparison. The change in 7 
exerts little effect upon the positive peak, but the 
negative excursion is delayed in its maximum by 
about 5 /xsec. The corresponding spectra are shown 
in figure 4. The spectral peak, / m , is lowered by 
about 1 kHz by the change in 7. 

The effect of variations of a and £ must be consid- 
ered simultaneously. A large number of spectra 
have been derived for strokes with 7=3 X10 4 sec -1 , 
^=3X10 8 m/sec; 1 , and F -=8X10 7 m/sec -1 , and 
various combinations of a and /3. Figure 5 shows 
f m as a function of a and k, where k=fi/a. It is seen 
that for large k, which corresponds to A^ery short rise 
times for the current surge, j m is determined almost 
entirely by a, that is, by the rate of decay of the 
current waveform. 

2.4. Return-Stroke Spectra Deduced From Current 
Measurements 

The various methods of measuring the current 
surges in return strokes have been described by 
Wagner and McCann [1950]. They include klydono- 
graphs, fusible wires, magnetic surge-crest ammeters, 
magnetic surge-front recorders, fulchronographs, and 
cathode-ray oscillographs. Care must be used in 



interpreting the data presented. For example, 
the magnetic link records give the maximum currents 
in flashes rather than in individual strokes. 

A review of data on return-stroke current wave- 
forms suggests that Berger 's [1961, 1962] measure- 
ments using a high-speed oscilloscope are of special 
merit. It appears that many previously reported 
values of rise times are in error because of instru- 
mental lag. This is especially true for subsequent 
strokes, which Berger has found to peak in roughly 1 
/xsec. For a first return stroke, transferring negative 
charge to earth, the observational data suggests the 
parameter values, a=2X10 4 sec -1 , 0=2 X10 5 sec -1 , 
and ^=30,000 A. Curve (a) of figure 6 shows 
d 2 M/dt 2 for such a stroke, with 7=3 X10 4 sec -1 , 
F =8X10 7 m/sec -1 , and ^=3X10 8 m/sec -1 . 

For subsequent strokes, the observations indicate 
that a equals 1.4 X 10 4 sec -1 , £is at 6 X 10 6 sec -1 , and i is 
10,000 A. For this case, k(P/a) is near 400, and so 
exerts little influence upon the VLF spectrum. 
Accordingly, (2.10) can be simplified by setting ft 
equal to 00. Curve (b) of figure 6 shows d 2 M/dt 2 
for this subsequent return stroke, with V and u as 
for the first return stroke and 7 equal to 0. The 
latter condition implies an infinite upward growth 
rate of the channel, but the error is not serious. By 
the time the length of the "model" channel surpasses 
I he actual one 4 , the current has already become very 
small. 

Flashes to earth bringing down positive charge 
are rare; they represent perhaps only 5 percent of 
the total number of discharges to ground. Berger 
[1962] reports that the strokes carrying positive 
charge to earth usually take about 30 /xsec to peak 
and that the maximum current is of the order of 
5,000 A. It seems that appropriate parameters 
would bei =7,000 A, a=SX10 3 sec -1 , and £=SX10 4 
sec -1 , while for lack of better data the values of V , 
7, and u used in the first stroke (negative) will be 
retained. The plot of d 2 M/dt 2 for this stroke is 
shown as curve (c) of figure 6. 

The spectra of the typical return strokes just 
described are shown in figure 7, where curves (a), 
(b), and (c) correspond to those of figure 6. Ft is 
seen that the peaks occur near 5, 3, and 3 kHz, 
respectively. 

Some writers have extended the Bruce-Golde 
model to include a third and even fourth exponential 
term, writing the expression for the return-stroke 
current as i t = i [a exp(—at) — b exp( — pt)-\-d exp 
(—8t)+e exp(-et)} [Barlow et al., 1954; Williams, 
1959]. Such expressions can be used to include the 
weak currents which sometimes persist for several 
milliseconds in lightning channels. Examination of 
oscillograms showing such currents reveals a tend- 
ency for them to be rather flat, with irregular surges 
(the M changes) distributed in a random fashion. 
They do not nt the exponential model, as a rule, and 
any attempt to find typical values of 8 and e does not 
appear to have much chance of success. The long 
continuing components are significant in the genera- 
tion of ELF radiation (slow tails) rather than VLF 
radiation. 
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Figure 1. Curves of d 2 M/dt 2 versus time at 100 km, showing effect of varying u. 
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Figure 2. Return-stroke amplitude spectra at 100 km, showing effect of varying u 



3. Observations of Return-Stroke 
Atmospherics 

The preceding section has shown how the source 
amplitude spectrum of theVLF atmospheric gener- 
ated during a return stroke may be derived by com- 
bining theory and the experimental measurements 
of current surges. Many researchers have also 
deduced amplitude spectra at the source by observing 
atmospherics. Since the atmospherics are usually 
recorded at a considerable distance from the origi- 
nating lightning discharge, it is necessary to apply 
an appreciable correction for the modifications to 
the spectrum introduced during propagation, if 
extrapolation back to source conditions is to be 
reliable. 



3.1. The Mean Amplitude Spectrum 

An examination has been made of the average 
amplitude spectrum as deduced from several re- 
searches. These include the work of the group at 
King's College, London, England [Chapman and 
Matthews, 1953; Edwards, 1956; and Steptoe, 1958] 
and other British papers [Horner, 1961; Croom, 
1961] and three American contributions [Watt and 
Maxwell, 1957; Taylor and Jean, 1959; Taylor, 
1963]. Altogether, 11 average curves were con- 
sidered. When these data are examined, after 
normalization to a distance of 100 km, it is found that 
the King's College results tend to be greater in 
amplitude than those of other workers. It is 
believed that this is due to the selection of the larger 
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Figure 3. Curves of d 2 M/dt 2 versus time at 100 km, showing effect of varying y. 
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Figuke 4. Return-stroke amplitude spectra at 100 km, showing effect of varying y. 
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a (10° sec ) 
Figure 5. Variation of f m with a and 0. 



atmospherics by triggering methods in an area of 
high background noise, and by an effective reliance 
upon a triggering system (the British CRDF net- 
work) for the location of the sources of the atmos- 
pherics. The observations of Horner [1961], in 
which the amplitudes at 6 kHz of all atmospherics 
originating in a local storm were recorded, are prob- 
ably by far the most reliable guide to the mean 
amplitude. 

Curve (a) on figure 8 shows the average derived 
from the experimental data considered. The plot is 
normalized in terms of E/E m and/// m where E is the 
field in /xVsec/m and/ the frequency. The frequency 
i m is that at which the peak value of field, E m , is 



attained. Also represented on figure 8 are the two 
equations 



and 



log(E/E m )=-\\ogf/j m \ 
log (E/E n ) = -(fl.8)\logflf m 



(3.1) 
(3.2) 



The value of j m lies between 3 and 7 kHz for the 
various groups of experimental results with 5 kHz 
being a good average; at 100 km, E m may be esti- 
mated as 50 juVsec/m, the greatest reliance for the 
amplitude information being placed on Horner's 
[1961] work. Scales corresponding to f m =5 kHz 
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Figure 6. Curves of d 2 M/dt 2 versus time at 100 km for typical return strokes. 



785 




=L 



20 



"l — I — rr 



v 



s 



N 



\ 



\ 



-(A) FIRST NEGATIVE STROKE 

-(B) SUBSEQUENT NEGATIVE STROKE 

-(C) POSITIVE STROKE 



\ 



\ 



\ 



\ 



\ 



J I L_X 



20 



FREQUENCY 



kilohertz 



Figxjee 7. Amplitude spectra at 100 km for typical return strokes. 



and E m =50 juVsec/m are included on figure 8. Over 
the limited frequency range of 2 to 20 kHz, (3.2) 
fits most of the experimental results reasonably well. 
Perhaps the worst agreement is with the data of 
Taylor and Jean [1959] for which the curve of E/E m 
against frequency is much broader than is indicated 
by (3.2); the Taylor-Jean measurements also give 
the highest value, 7 kHz, for j m . This latter result 
is perhaps not unexpected, since the Taylor- Jean 
data included many mountain storms for which there 
is reason to believe the channels are short and y is 
large; hence (sec. 2.3) j m is high. It is interesting 
to note that Taylor's [1963] most recent experimental 
data, obtained in the Texas-Oklahoma area and 
therefore uncontaminated by mountain storms, give 



f m for return-stroke atmospherics as averaging 5 
kHz. This is in excellent agreement with the general 
conclusions above. 

3.2. Amplitude Variations 

Variations in amplitude, generated at the source, 
may be observed upon a broadband or narrowband 
basis. If the source spectra of all atmospherics were 
the same, which is obviously an unrealistic postula- 
tion, identical results for amplitude distribution 
would be obtained by the two approaches. Valuable 
results regarding source amplitude variations are 
best obtained under two sets of circumstances. 
First, a close storm may be studied; in this instance 
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Figure 8. Normalization of average spectrum for return-stroke atmospherics. 



propagation effects are negligible. Second, atmos- 
pherics from an isolated distant storm may be 
recorded over a period during which propagation 
conditions are stable. 

An examination of several samples of data shows 
that in all cases the amplitude variability is well 
represented by a log-normal distribution. If E^ed is 
the median amplitude and the size 2? of the individual 
atmospherics is expressed in decibels relative to 
^med, then P, the percentage exceeding E p , is given 

by 

100 f 00 / ff-£ med y^ 

fik^-K-^)**' (8 - 8) 



p-- 



T V27 



where a is the standard deviation in decibels. On 
probability paper, the relation between P and E is a 
straight line with a slope of a. The log-normal 
distribution is, of course, that expected where a 
number of contributing factors of independent vari- 
ability enter into a problem. 

A typical sample of observations is plotted on 
fig. 9. These are for the peak-to-peak amplitude 
of radiation fields from a storm at a distance of 
60 km [Pierce, 1950]. A straight line, corresponding 
to the log-normal distribution, fits the results quite 
well, the slope being 8.6 dB. Other information 
that has been examined also gives log-normal varia- 
tions. Among these data are the 10-kHz observa- 



tions of Chapman [1955] ((7=6.7 dB) and the 10.5- 
kHz measurements of Hefiey, Doherty, and Linfield 

[1959] 0=4.0 and 6.3 dB). Although the radiated 
atmospheric is proportional to d 2 M/dt 2 , where M 
is the electric moment and t the time, and should 
therefore not necessarily have the same amplitude 
distribution as variations in current or M, it is 
perhaps worth noting that the peak current magni- 
tudes listed by Bruce and Golde [1941] fit a log- 
normal distribution with a =8.5 (IB. Furthermore, 
two samples of data given by Pierce [1950] for the 
change in M associated with a return stroke also 
follow the log-normal distribution, with standard 
deviations of 6.5 dB and 9.5 dB. All these results 
are in reasonable agreement with the conclusion 
of Horner [1960] that the amplitude variations of 
return-stroke atmospherics at the source obey a 
log-normal law with a standard deviation of some 
7dB. 

3.3. Spectral Variations 

It has been tacitly suggested in the preceding two 
sections that the shape of the amplitude spectrum 
is defined by (3.2), that spectral fluctuations are 
taken account of by changes in f m , and that ampli- 
tude variations apply similarly at all frequencies 
with an overall multiplying factor B that follows 
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Figure 9. Amplitude distribution for peak-to-peak radiation fields at 60 km. 



a log-normal distribution, 
modified to the form 



Equation (3.2) can be 



E=BE m antilog{-(0.8)|log/// w 



(3.4) 



It follows that if amplitudes are being compared at 
two frequencies, fi and f 2 , both exceeding f m , then 
the ratio E n /Ef 2 is independent of both f m and B, 
and apparently invariant for two given frequencies. 
If, however, a sample of atmospherics is examined in 
this manner, the ratios are not found to be constant 
but to follow a distribution that again approximates 
a log-normal law. 

Figure 10 shows an illustrative sample. The 
data are those of Chapman [1955], the ratios of 
amplitudes at 18 and 10 kHz being recorded for 
atmospherics propagated over the same distance 
with similar ground and ionospheric conditions. 



The log-normal distribution is followed, with a =5.2 
dB. The results of Bowe [1951] also fit similar 
distributions, although the condition of j x and f 2 
being greater than j m is hardly satisfied for Bowe's 
samples. At 4.25 and 7.5 kHz, (7=6.2 dB; at 3.5 
and 7.5 kHz, two samples yield o-=4.8 dB and <7=5.2 
dB. 

The apparent discrepancy between the observed 
fluctuations in the ratio of amplitudes at two 
frequencies and the constancy predicted by (3.2) 
is essentially due to the deficiencies of (3.2). Al- 
though it is perhaps the best simple formulation 
of the mean spectrum for a group of atmospherics, 
and the differences in the mean between groups 
may be represented by adjusting f m , it is quite 
incorrect to consider that (3.2) should fit all indi- 
vidual atmospherics. 
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Figure 10. Amplitude distribution of ratio 18 to 10 kHz. 



4. Discussion and Summary 

In section 2, calculations on a semi theoretical 
basis were made and the amplitude spectra at 100 
km deduced from the current surge information. 
In section 3 a comparable spectrum was derived 
from observations of atmospherics. For comparison, 



figure 11 shows the three theoretical curves ((a), (b), 
and (c)) obtained for a first stroke, a subsequent 
stroke, and a positive stroke, respectively; also given 
are the average obtained from the experimental 
atmospherics observations (curve (d)) and (3.2) 
(curve (e)) which has been shown in section 3 to be a 
good fit over the frequency range of 2 to 20 kHz to 
the mean of most sets of experimental results and 
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Figure 11. Summary of amplitude spectra at 100 km. 



which, with E m =50 /xV sec/m and / w =5 kHz, is an 
excellent representation of the overall average. 

The agreement in figure 11 between curves (a), 
(b), (d), and (e) is very good. The observation of 
atmospherics may be expected to sample both first 
and subsequent strokes; accordingly (d) should be 
intermediate between (a) and (b) both in amplitude 
and in the value of j m . This is approximately what 
is found. Experimental data on positive strokes 
are not readily available, but by analogy with the 
negative strokes, curve (c) may be expected to be a 
good representation. 

The consistency found between the theoretical 
approach and the spectrum obtained by extrapolation 
from atmospherics, with the amplitude adjusted 



ignoring the King's College results, is evidence for 
the belief advanced in section 3 that this school of 
workers overestimates the size of the average at- 
mospheric. It is shown in appendix 1 how such 
overestimation can result from the selectivity 
intrinsic in the use of triggering methods; these 
methods can also lead to incorrect deductions re- 
garding propagation even when only one waveguide 
mode of propagation is concerned. 

The range of j m from the experimental observations 
was between 3 and 7 kHz. This is entirely in agree- 
ment with the variation shown in section 2 to result 
from changes in the parameters involved in the 
theoretical approach. An apparent paradox is that 
j m is less for subsequent strokes than for a first 
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Figure 12. Variation of amplitude with frequency for f>f m . 



return stroke, whereas the current surge is more 
sudden in the former case. The explanation of the 
paradox is, of course, that the length of channel as 
well as the current is involved in determining the 
radiated spectrum. In the limit, an infinitely rapid 
current pulse would not radiate any atmospheric, 
since the channel length (effective radiating antenna 
length) would be zero. 

Amplitude variations obey a log-normal distribu- 
tion with a standard deviation of between 4 and 10 
dB; 7 dB is an average value. As regards spectral 
variations, beside the dissimilarities between groups, 
represented by an adjustment of f m , considerable 
differences exist between separate atmospherics. 
Notably, the ratio of the amplitudes at two selected 
frequencies follows a log-normal distribution, with 



o-=5 dB. This kind of variation can be accounted 
for on the theoretical approach by changing the 
parameters and introducing subsidiary current 
surges. However, the data presently available are 
far too restricted to justify a treatment of these 
secondary influences with any degree of confidence. 
The variation of amplitude with frequency for 
frequencies exceeding/™ has been a question given some 
attention in the past. Much effort has been devoted 
to fitting the variation to a law of the form f~ c 
where c is a constant. It is now becoming evident 
that a particular value of c is only applicable over a 
limited range in frequency and that as this range 
increases so does c. For example, the experimental 
data of figure 8 give c = 0.S between 5 and 20 kHz; 
the results of Taylor [1963] show c=1.0 for the 



791 



frequency range of 10 to 100 kHz and somewhat 
less from 5 to 10 kHz; while Watt and Maxwell 
[1957] suggest c = 2 between 40 and 80 kHz. The 
theoretical treatment also implies c increasing with 
frequency. This is illustrated in figure 12 in which 
the spectrum of figure 11 for a typical first stroke is 
plotted on a log-log basis over the range from 1 to 
100 kHz; for comparative purposes the lines cor- 
responding to c = 0.5, 1.0, 1.5, and 2.0, and passing 
through the point (E m / m )are shown. 

In conclusion it is perhaps appropriate to stress 
the dangers of extrapolating laws applying to the 
VLF range to higher frequency bands, particularly 
if the general radio noise from a thunderstorm and 
not that from particular types of atmospherics is 
being considered. Even for the return stroke the 
theoretical development of section 2 is very idealized: 
no account is taken, for example, of branching, a 
feature which may well be responsible for some of 
the spectral irregularities experimentally observed 
at LF [Taylor, 1963]. If the total noise generated 
by lightning at VLF and LF is being considered 
then obviously the disturbances due to the leader 
and to K pulses, and discussed in the companion 
paper [Arnold and Pierce, 1964], should be included 
in addition to return-stroke atmospherics. Even 
so the description of the noise may not be complete: 
small sparks not considered in either of the two 
present papers are undoubtedly the dominant noise 
generators at HF [Pierce, 1962], and may also have 
some significance at lower frequencies. The ulti- 
mate limits in thunderstorm noise generators are 
the minute sparks associated with corona and other 
discharges from water drops [Pierce, 1962; Sartor, 
1963]; this noise source may bear little relation to 
the actual occurrence of a lightning flash. 

5. Appendix 1 

Very-low-frequency atmospherics produced by 
leader pulses, by K changes, and by return strokes 
in a given thunderstorm are necessarily recorded 
against a background of circuit noise, manmade 
interference, and interference from other storms 
occurring at the same time. To reduce such inter- 
ference many recording devices utilize triggering 
mechanisms, so that only those atmospherics exceed- 
ing a predetermined threshold amplitude are recorded. 
This means that the mean intensity of the recorded 
strokes from distant storms, considered as a function 
of range, does not decrease as rapidly as the attenua- 
tion function, taken by itself, would indicate. 

Section 3 indicates that the source amplitudes of 
individual atmospherics from a nearby storm follow 
a log-normal distribution. That is 
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where E is the amplitude of an individual atmos- 
pheric expressed on a logarithmic scale; E mea is the 
median amplitude; Q{E) is the probability of an 



individual atmospheric exceeding amplitude E; 
and a is the standard deviation for the sample. 
Experimentally, a was found to be between 4 and 
10 dB. 

Let the threshold of the detection device be E . 
The mean amplitude of the recorded atmospherics, 
measured with respect to the threshold, is 
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This function is shown versus (E —E mea ) in 
figure 13, plotted in units of a. For cases when 
practically all the atmospherics are recorded 
(E^^^med), raising the threshold by a given number 
of decibels reduces E by almost the same amount. 
As E becomes larger than E med , however, the change 
in E for a given change in E (or E mea ) becomes quite 
small. 

It is necessary in analyzing the amplitude distri- 
bution of a sample of atmospherics to take into 
account the relationship between E and E mea ; E 
may be set to limit the number of strokes recorded 
or to avoid triggering by random noise. Circuit 
noise in the receiving equipment proves to be 
negligible, but manmade noise picked up by the 
antenna is significant in some areas, especially near 
large cities. A typical value for E in the centra] 
London area, where much work on atmospherics 
has been done, is 2 iiV sec/m at 10 kHz. Comparing 
this with the field produced by a typical return 
stroke 100 km distant of, say, 50 fxV sec/m shows 
that even under such unfavorable conditions the 
signal-to-noise ratio for the median stroke is around 
25 to 30 dB for a nearby storm. Taking a, the 
standard deviation of E, as 7 dB, it is seen that the 
proportion of strokes masked by noise is completely 
negligible, being of the order of 10~ 5 . 

Severe manmade interference can distort observa- 
tions on distant storms. The effect of attenuation 
is to shift E mea down, with a remaining the same. 
The attenuation in decibels at 10 kHz is given 
approximately by [10 log d+(2.5)lO~*Xd\ where d 
is the distance in km. At 3 000 km, an atmospheric 
is reduced in amplitude 22 dB below its amplitude 
at 100 km. This implies a signal-to-noise ratio 
near one for a median stroke recorded at 3000 km 
in a noisy area, with half of the strokes not recorded. 
This situation corresponds to the region near 
(E —E mea )=0 on figure 13; some distortion in 
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Figure 13. E as a function of threshold setting. 



experimental determination of E is to be expected 
in such a case. 

In practice, the tendency in recording atmospherics 
is to set E high enough so that only well-defined 
waveforms are recorded. Since these must stand 
out, not only against any manmade noise that may 
be present but also against a persistent background 
of weaker or more distant atmospherics, only the 
strongest atmospherics are recorded. For example, 



some studies in Great Britain have been limited to 
the atmospherics strong enough to be located by the 
Meteorological Sferics Network [Ockenden, 1947]. 
Such records can yield considerable information, but 
the limitations implicit in triggering methods must 
be borne in mind. 

The use of such sampling techniques necessarily 
involves working in a region toward the right side 
of figure 13, that is, with (E —E med ) being equal to 
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2(r to 3 cr. In this case, there is little variation in E with 
changes in E or E med . Therefore, E is a very 
slowly varying function of d, and sampling errors 
will largely determine which, of two distant storms 
yields a greater value of E. This is a plausible 
explanation for the apparently negative attenuation 
coefficients that have been reported by some 
researchers. 



We are grateful for the use of data unpublished 
in the open literature but contained in the theses 
of Dr. B. J. Steptoe, Dr. J. Chapman, Dr. D. L. 
Croom, and Dr. J. C. Williams. 

The research described in this paper summarizes 
some of the work discussed in greater detail in a 
report by Pierce, Arnold, and Dennis [1962]. This 
report was prepared under U.S. Air Force Contract 
AF33(657)-7009; it is a pleasure to acknowledge the 
action of the Air Force authorities in permitting 
publication. 
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